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The heavy-quark QCD sum rules are applied to a model of charmonium based upon the 
gauge/gravity duality. We find that there is strong agreement between the moments of the po- 
larization function calculated from the holographic model and the experimental data suggesting 
that the model is consistent with the heavy-quark QCD sum rules at zero temperature. 
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I. INTRODUCTION 

A new model of charmonium in the context of the 
gauge/gravity duality was constructed at zero temper- 
ature using spectral data so that the model exactly re- 
produced the masses and decay constants of the ground 
state, J/ip, and the first excited state, ip' [l|. This model 
was then considered at finite temperature, where the dis- 
sociation temperature of charmonium was investigated. 
In this paper, we will examine the zero temperature 
model in the context of heavy-quark QCD sum rules 0- 
4]. As a reference point, the results will be compared 
with a similar analysis of a different holographic model 
of charmonium [B|, Q . 

QCD sum rules were first developed by Shifman, Vain- 
shtein, and Zakharov (SVZ) @, HI]. They provide a 
systematic way to compare theory and experiment and 
have been used in a wide variety of systems. QCD sum 
rules are based upon the two-point vector-vector corre- 
lation function of the heavy-quark (charm) current; at 
zero temperature this correlation function can be ex- 
pressed in terms of one function, the polarization func- 
tion, lK c )(Q 2 ), 



Ae" 1 <Jn(x)JM) = (9/»9* - <?V) n (c) (Q 2 ), (1) 



where the charm current is defined as J M = cj^c and 



Q 



-q . The charm polarization function can further 



be expressed in terms of a dispersion relation 

_ _^_ n ( c )/ 2N = 1 f Rc(s)ds 
dQ 2 W) 12n 2 Q 2 J (s + Q 2 ) 2 '' 



(2) 



where Q c = 2/3 is the charm quark electric charge, and 
R c (s) is the imaginary part of the polarization function 
and is related to the cross section a c as 



R c 
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(3) 



QCD sum rules states that the cross section, and 
thereby the polarization, can be formulated either from 
experimental data or theoretically from an operator 
product expansion (OPE). Since both of these should 



describe the same physics, the polarization function one 
obtains from either method should be the same. This 
leads to the construction of equations between the po- 
larization function calculated using either the OPEs or 
the experimental data. These equations can be solved for 
various parameters of the OPEs. 

As pointed out by SVZ, the OPE side of the sum rule 
becomes tractable when one considers a regime where 
asymptotic freedom allows one to use perturbative QCD. 
In general, this is valid for energies which satisfy Q 2 + 
^m 2 q 3> Aq CD . For light quarks, this is only satisfied 
when Q 2 Aqcd- However, in the case of heavy quarks, 
such as charm, the relation holds even when Q 2 = 0. 
SVZ used this fact to construct QCD sum rules which 
arc unique to heavy-quark systems 1 . In the region of 
small Q 2 , it is useful for comparisons to introduce the 
moments of the polarization function defined as 
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n< c >(Q 2 ) 



(4) 



Q 2 =0 



One can then relate each moment between the OPE side 
and the experimental data rather than the entire polar- 
ization function. In the orig inal work by SVZ using 
this method of comparison, they were able to calculate 
the heavy-quark mass and show the need for nonper- 
turbative corrections to have better agreement with the 
data. 

Recently, the development of gauge/gravity dualities 
based upon the AdS/CFT correspondence [8l4l0l| have 
led to holographic models which provide a new method 
of calculating the polarization function. Therefore it is 
natural to apply QCD sum rules as a way to understand, 
constrain, and evaluate possible holographic models. The 
OPE expansion has already been used for light quark 
holographic systems in order to show that the AdS met- 
ric is sufficient to reproduce the leading large Q 2 behav- 
ior [ll|, [l2| , and agreement with the OPE expansion has 
been a guiding principle in the construction of ma ny l igh t 
quark holographic models among which includes [l3l - [r7| . 
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1 Though in this note, we focus on the heavy quark sum rules at 
n2 —■ 0, heavy quark sum rules have been used at Q 2 > based 



upon the initial work of Ref. Q 
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Yet this is the first time that these techniques will be 
applied to a holographic model with heavy quarks. 

In the next section, the gauge/gravity dual models of 
charmonium 0, [j| are introduced. In Sec. IIII1 the mo- 
ments of the polarization function will be calculated from 
perturbative QCD as well as the holographic model 
This will be followed in Sec. IIVI by the results of the 
comparison of the moments between QCD and the holo- 
graphic models, where the holographic models will take 
the role of the OPE side in the sum rules. Finally, we 
will conclude, in Sec. E] with a discussion of the results. 



II. CHARMONIUM MODEL FROM 
GAUGE/GRAVITY DUALITY 

Key features relevant for the current discussions of 
holographic models of charmonium will be presented 
here. A more detailed exposition of these models can be 
found in Refs. 0] and Q. The model of Ref. 0] will be re- 
ferred to as the "shift and dip" model while the model of 
Ref. @ will be referred to as the "rescaled p" model. The 
basic construction of the two models is similar, but the 
latter does not reproduce as many phenomenologically 
relevant features as the shift and dip model. The differ- 
ences between the two models will be pointed out, and 
the two models will each be compared with the heavy- 
quark QCD sum rules in the subsequent sections. 

In the spirit of the holographic approach, it is assumed 
that the generating functional of the heavy-quark vector 
current J M can be represented by the effective action ob- 
tained by integrating over a bulk 5D gauge field Vm (dual 
to the current) at a given fixed boundary value (equal to 
the source of the current). The action for the 5D gauge 
field is given by 



1 

^1 



S = —T2 I d & x^e-"V M NV MN , 



(5) 



where g 2 is the 5D gauge coupling and Vmn — 9mVm — 
9nVm- The two-point current correlator is given by the 
linear response of the field Vm to an infinitesimal pertur- 
bation of its boundary condition. 

The conformally fiat representation for the 5D back- 
ground metric gMN with 4D Lorentz isometry is chosen: 

ds 2 = g MN dx M dx N = e 2A{z) [r) fiU dx>*dx v - dz 2 ] , (6) 

where rj^ = diag(l, — 1, — 1, — 1) is the Minkowski met- 
ric tensor. The effect of confinement is represented by 
the nontrivial background profile of the scalar field 3> in 
Eq. (|S|) in the same way as it is done in the soft-wall 
model with dilaton background in Ref. (l8| . 

Following the rules of the holographic correspondence, 
the generating functional for correlation functions of the 
heavy-quark current can be calculated by evaluating the 
action at its extremum for given boundary conditions. 
The extremum is given by the solution of the equations 
of motion, which in V5 = gauge read 

Z [e B ^d z V} + q 2 e B( - z V = , (7) 



where V is any of the three components V± of V^(z,q) 
transverse to 4-vcctor (q 1 - = 0) and 



B = A- $. 



(8) 



Discrete values of q 2 — m 2 , for which Eq. ((7} possesses 
a normalizable solution V — Vk{z) satisfying the bound- 
ary condition 2= o = 0, correspond to the masses of 
the charmonium states, k = 1,2,... = J/i/j,ip' , . . . . We 
normalize such solutions as 



dze B{z) v k {z) 2 = 1 



(9) 



The current-current correlator can be calculated ac- 
cording to the well-known prescription of Ref. pjl, [2(| , 
by 



G R {q) 
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2 e B V(z,q) 



(10) 



where e — > is an ultraviolet regulator and V{z, q) is 
the non-normalizable solution of Eq. (J7J with boundary 
conditions: 



= 1 

Z— >oo 



V{e,q) 

V(z,q)^^0 
The polarization function can then be determined by 



(11) 



Tl^(Q 2 )=--^e B V'(z,q) 



(12) 



where Q 2 = —q 2 . The ultraviolet behavior of the 
current-current correlator is required to be conformal, 
viz. Gn(q) ~ q 2 \og{—q 2 ) as q 2 — > —00, which translates 
into 



(13) 



and matches that of QCD, which fixes [111. Il2t |21| 

eg = Un 2 /N c . (14) 
By performing a Liouville transformation, 

we can bring Eq. ((TJ) to the canonical Schrodinger-like 
form, 



-d 2 V/dz 2 + U(z)iS> = qH 
with the holographic potential given by 



U{z) 



B"{z) (B\z) 



(16) 



(17) 



At this point the two holographic models differ. The 
rescaled p model chooses the metric warp factor to be 
A(z) — — log(z) and the dilaton profile to be $(2) = 
a 2 z 2 . This is precisely the form of the standard soft-wall 
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model [18]. This leads to a model of one parameter, viz. 
a, which is determined by fixing the mass of the J/ip. 
Therefore the rescaled p model correctly reproduces the 
J/ip mass (by construction), but the J/ip decay constant 
and the ip' mass and decay constant are all off by nearly 
20%. 

On the other hand, in the spirit of the bottom-up ap- 
proach, the shift and dip model chooses the function B 
so as to satisfy the spectroscopic data associated with 
J/ip and ip 1 . It is assumed that such a background arises 
dynamically, but no attempt to model the corresponding 
dynamics is made. To this end, a holographic potential 
U(z) is chosen: 



+ ((a 2 z) 2 + c 2 ) 0(z - z d ) - a5(z - z d ). 



(18) 



There are four parameters in this potential, which will 
be used to fit four experimental data points: the masses 
and the decay constants of J/ip and ip'. This potential 
may seem rather unusual, but it was chosen to reproduce 
the correct expected behavior at both small and large 
z as well as adding features, viz. the c 2 "shift" term 
and the delta function "dip" term, which facilitated the 
matching of the parameters to the physical constraints. 
The matched parameters were found to be 



a = 0.970 GeV, c = 2.781 GeV, 
a = 1.876 GeV, z~ x = 2.211 GeV. 



(19) 



The potential in Eq. (|18p can be, certainly, improved 
by applying further constraints. One goal of this paper 
is to determine how well this potential holds up to the 
additional constraints imposed by the heavy-quark QCD 
sum rules. 

Having established U(z), the polarization function can 
be found by first solving Eq. (ITBl for the function ^(z). 
Equation (jT7)) can be solved for the function B(z). There- 
fore the vector field can be calculated from B(z) and 
$f(z) using Eq. (|T5j) . Finally the polarization is calcu- 
lated from Eq. (|T2"j) while its moments can be generated 
from Eq. (gj). 



A. OPE expansion 

To calculate the moments for the QCD OPE expan- 
sion, we will concentrate on the leading unit operator 
term. Nonperturbative corrections, due primarily to the 
gluon condensate, can be included, but we will later show 
that this is unnecessary for our analysis here. This term 
can be calculated using perturbative QCD because of the 
heavy charm quark. At the one-loop level, the imaginary 
part of the polarization function is given by 



Rc(s) = -v(3 



v 2 ) 6>(s-4m 2 ), 



(20) 



where the heavy-quark velocity is given by v = (1 — 
4m 2 /s) 1 / 2 , and m c is the charm quark mass. This ex- 
pression can be used in the dispersion relation of Eq. (j5J) 
to calculate the polarization function, or it can be used 
in Eq. ©to find it moments. This leads to the familiar 
result of [3j for the moments, 



1 



3 2"(n +!)(»- 1)! 
4^2 (2n + 3)!! (4m 2 f 



(21) 



where superscript (0) refers to the fact that this is the 
leading order calculation. Note that charm quark mass 
is a free parameter of the OPE. 

The leading a s correction to the moments can also 
be calculated. With the inclusion of this correction the 
moments are given by 



(22) 



where 



/ VtF r(n + 3/2) 1 - (3re + 3)~ ] 
V~3 T(n + 1) l-(2n + 3)- ] 
3 4nlog2 

7T 

3 r(n + 3/2)l-2(3n + 6)~ 

3Vtt 4tt T{n + 2) 1 - (2n + 3)" 
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(23) 



B. Shift and dip model 



III. CALCULATION OF MOMENTS 

In this section, the moments will be calculated directly 
from the polarization function for the QCD OPE expan- 
sion and the shift and dip holographic model. The cal- 
culation of the polarization function using the rescaled 
p model is presented in Appendix The method to 
calculate the same moments from spectroscopic data will 
also be explored. A more general calculation of the mo- 
ments for an arbitrary holographic model can be found 
in Appendix [Bl 



In this section, the procedure discussed at the end of 
Sec. im will be used to determine the polarization function 
for the shift and dip model. The moments will not be 
explicitly calculated, but can be found from Eq. (j4j using 
the polarization function. The results will be expressed in 
terms of the parameters a, c, a, and z d whose numerical 
values can be found in Eq. (Til?]) . First, one needs to 
calculate B(z); this is done by solving Eq. (fT7| with the 
holographic potential giving by Eq. (fT8"|) , This results in 



,B{z)/2 



7t + Cl 



c 2 2" 



z < z d 
H V1 (az) z > z d 



(24) 
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where H v {x) is the Hermite polynomial, v\ = — | — js, Eq. (fl6]) and use Eq. (|T5|) . This results in 
and ci and ci are integration constants given by, 



(l + 2az d + 2a 2 zj) H V1 (az d ) - 2az d H V2 (az d ) 
(3 - 2az d - 2a 2 zj) H Vl (az d ) + 2az d U v , 2 (az d ) ' 



c 3 (Q 2 ) qzJ 1 (qz)-% qz Yj(gz) 



>> •>•>-■ '>" 2 ^)ff,>z d ) + 2 a z d I^ 2 (az d )' 4) = J (26) 

JZdH Vi (az d ) v 



(25) 

where z^i is as before and 1^2 = ^1 + 1. In order to de- where c\ and C2 are as before and c^{Q 2 ) and Ci(Q 2 ) are 
termine the bulk vector field, V{z, q), one needs to solve given by 

I 



,q 2 \ — Hi ^ 2zd( i Y ^ z d> H ^M z d) + Y i(qzd) ((1 + 2az d + 2a 2 zj) H viq (az d ) - 2az d H U2q (az d )) \ 
2 y -2z d qJ (qz d )H„ lq (az d ) + Ji(qz d ) ((1 + 2az d + 2a 2 z 2 d ) H Vlq (az d ) - 2az d H U2q (az d )) J 

(n 2^ _ H Vl {az d ) c 3 (Q 2 ) gzdj^qzd) ~ \qz d Y x ((iz d ) 
tW j ~ H vi (az d ) 1 + ci 



with q 2 = -Q 2 , and V\ q = V\ + and v 2q = v Xq + 1. 
From this formula for the vector field, the polarization 
function can be found using Eq. (TT2"1) . 

n W( ^ ) = ^( TC+1 „ g (f)- ca (g V ^). 

(28) 

One can show that, for small Q 2 , there are terms in 
c i{Q 2 ) which cancel both the c\ term and the log(Q) 
term. This expression also has the standard log(e) di- 
vergence from the UV cutoff. This can be removed by 
renormalization such that I[( c )(0) = 0. The moments can 
be generated from this expression using Eq. (Q| , however 
a simple analytic expression is not possible and therefore 
will not be presented here. 



C. Spectroscopic data 

The moments of the polarization function can also be 
determined directly from the spectroscopic data. If one 
assumes that the cross section can be constructed as a 
series of delta functions as is done in Refs. [H, HJ|, one for 
each charmonium state, R c can be written as 

16tt 2 f 2 (29) 

6 k k 

where the sum is over all charmonium states with k la- 
belling each state, T ee is the electronic width, and / is 
the decay constant. The moments can then be calculated 



by using this formula in Eq. (f?]), 

Since the mass of the excited states tends to grow faster 
than its decay constant, the higher excited charmonium 
states are suppressed for the higher moments. When 
considering the experimental spectrum, the spectrum ex- 
hibits a continuum of states above the DD threshold. 
SVZ approximated this continuum by choosing R c (s) to 
be 

R c (s) = ^(.s-16GeV 2 ), (31) 

above the threshold. The contribution to the polarization 
function and the moments can as before be found by 
using Eq. (|31~T) in Eqs. ([2]) and (HI), respectively. 

The necessary spectroscopic data, viz. the masses, rrik 
and decay constants, can either come from the physi- 
cal experimental results 22] , or it can be calculated from 
the spectrum of either holographic model. Unlike the ex- 
perimental data, the holographic models contain only an 
infinite tower of zero width excited charmonium states 
with no continuum. Therefore, for good agreement be- 
tween the moments of the holographic models and the 
experimental spectrum, one would hope that this tower 
of states can approximate the continuum of QCD. 

The equivalence between the two methods to calcu- 
late the moments, i.e. from the polarization function 
or from the spectrum, is due to equivalent ways of ex- 
pressing the polarization function. One tends to use the 
method which is easier to calculate with. In this paper, 
discussions will highlight the method using the spectro- 
scopic data since it is more natural to interpret the results 
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in the language of this method, though both techniques 
yield the same numerical results. 



IV. RESULTS 

To assess the consistency of the holographic model 
of charmonium with the heavy-quark QCD sum rules, 
we will examine the moments of the polarization func- 
tion qualitatively and quantitatively. The moments cal- 
culated from the holographic model will be compared 
to those calculated from the experimental spectroscopic 
data and from perturbative QCD. In order to provide a 
baseline for the results of the shift and dip model, the 
results using the rescaled p model will also be presented. 

To begin, we will use the techniques of SVZ to an- 
alyze the moments of the polarization function. The 
polarization function's moments will be calculated for 
each holographic model using each model's spectral data 
as detailed in Sec. IIII CI For both holographic models, 
a sufficient number of excited states were considered 2 . 
The moments of the polarization function are also calcu- 
lated from the first two experimentally known charmo- 
nium states and a contribution for the continuum as done 
in Q. 

The ratio of the moments of the polarization function 
between each holographic model and the experimental 
data are constructed and presented in Table U with su- 
perscript A referring to the shift and dip model and su- 
perscript B referring to the rescaled p model. The results 
indicate that the shift and dip model is consistent with 
QCD at worst at the 20% level with even better agree- 
ment for the higher moments. This is quite reasonable 
considering that there should be 1/N C corrections to the 
holographic model. 



Moment n 




•^„ S /^„ Expt 


1 


0.79 


0.54 


2 


0.91 


0.57 


3 


0.96 


0.59 


4 


0.98 


0.60 


5 


0.99 


0.61 



TABLE I: Ratio of the polarization function moments cal- 
culated using either holographic model to the moment cal- 
culated using QCD. The label A refers to the shift and dip 
model while the label B refers to the rescaled p model. 

The better agreement between the shift and dip model 
and QCD for the higher moments can be understood by 
the fact that the higher excited states are suppressed for 
these moments as discussed in Ref . [3|, |4( . Therefore the 



2 For the shift and dip model, 400 states were used, while for 
the rescaled p model, one can analytically calculate for infinite 
number of states. 



largest contribution for both the holographic model and 
QCD comes from the J/ip and tp' states which are identi- 
cal in both cases by construction. For example, the J/ip 
and ip 1 contribute only 75% to the first moment in the 
shift and dip model while this climbs to 99% for the fifth 
moment. Therefore the larger discrepancy seen in the 
first moment is an indication that the tower of excited 
states in the holographic model only marginally approx- 
imates the QCD continuum. 

The rescaled p model also has improved agreement 
with QCD for the higher moments because of the increase 
in the contribution of the J/ip to these moments. How- 
ever, the differences seen between the shift and dip model 
and the rescaled p model in Table Q] can be attributed to 
fact that the rescaled p model does not accurately repro- 
duce the ip' state. 

Another useful way to compare the holographic model 
to the QCD OPE and the experimental data is to look 
at the ratio of consecutive moments, defined as 



Figure Q] plots the results for r n from perturbative QCD, 
the experimental data, and both holographic models. 
The moments of the polarization function for the exper- 
imental data and the holographic models are calculated 
as described above. The polarization function's moments 
from perturbative QCD are found from Eq. (|23[) with 
ttlq = 1.21GeV. The charm quark mass was found by 
equating the first moment calculated from the QCD OPE 
and the first moment calculated from the experimental 
data and solving for the heavy-quark mass. The points 
at n = 1 correspond to the value of the first moment 
multiplied by the factor (47r/3) 2 , as is done by SVZ g|. 
Again the discrepancy seen between these values is in- 
dicative that the tower of excited states of the shift and 
dip model does not reproduce the experimental contin- 
uum well. One can see that for the highest moments the 
holographic models agree significantly well with the ex- 
perimental data. Again this is not that surprising since 
the highest moments have a larger contribution from the 
J/ip state. The discrepancy seen between the shift and 
dip model and the rescaled p model at the highest mo- 
ments is because the rescaled p model only reproduces 
the J/ip mass and not the decay constant. Furthermore, 
the perturbative QCD calculation with only the leading 
a s corrections, does not agree very well with the experi- 
mental data for the higher moments. SVZ showed Q that 
this discrepancy can be resolved by including the effects 
of the gluon condensate. The good agreement between 
the holographic models would imply that they correctly 
capture the nonperturbative physics of the charm system. 

The quantitative results from the moments can be 
summarized graphically by plotting the polarization 
function directly. Figure [2] depicts the polarization func- 
tion as a function of Q 2 for small Q 2 calculated using the 
experimental spectroscopic data (solid curve), the shift 
and dip model (dashed curve), and the rescaled p model 
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FIG. 1: Ratio of consecutive moments calculated from pertur- 
bative QCD (squares), spectroscopic data (circles), the shift 
and dip model (diamonds), and the rescaled p model (trian- 
gles). 



(dotted curve). Each function has been renormalized so 
that n( c )(<3 2 = 0) = 0. The curve from the experimen- 
tal data is calculated from Eq. ([2]) using Eqs. (|29|) and 
PTj) . The shift and dip curve is a plot of Eq. (gS]) with 
the parameters chosen from Eq. (fT9"| while the rescaled 
p curve plots Eq. (| A3[) ■ Note, though not shown, the 
polarization function determined from the perturbtative 
QCD calculation would lie directly on top of the exper- 
imental curve. From the figure, one can clearly see that 
the charmonium model with a shift and a dip has better 
agreement with experimental data as compared with the 
rescaled p model. The polarization function calculated 
from the shift and the dip model is in modest agreement 
with the QCD calculation (within 21%) for Q 2 < lGeV 2 , 
but remains within 26% agreement to Q 2 = lOGeV 2 . 
This uncertainty is associated with the discrepancy of 
the first moment seen in Table |TJ 




Q 2 (GeV 2 ) 



FIG. 2: Polarization function Il' c ' at small Q 2 calculated from 
the experimental data (solid line), the shift and dip model 
(dashed line), and the rescaled p model (dotted line). All 
functions have been renormalized so IT ' (0) = 0. 



V. DISCUSSION 

In this paper, we have considered a holographic model 
of charmonium and analyzed it in the context of heavy- 
quark QCD sum rules. We have illustrated how to cal- 
culate the moments of the polarization function for this 
model. The moments have then been used to compare 
the holographic models to the results obtained from the 
QCD OPE expansion and the experimental spectroscopic 
data. 

The results have shown that the holographic model of 
charmonium presented in Ref. [l| agrees very well, more 
than 80% for the first moment and reaching above 90% 
for the higher moments, with heavy-quark QCD sum 
rules. This was shown via graphical means and by di- 
rectly examining the moments of the polarization func- 
tion. 

This has been the first time that heavy-quark QCD 
sum rules have been applied to a holographic model. 
It is not obvious a priori that for the heavy-quark sec- 
tor the agreement between the holographic model and 
QCD for Q 2 < translates into agreement between the 
holographic model and QCD for Q 2 > 0. Nevertheless, 
within reasonable uncertainties, this is exactly what we 
have demonstrated here. The results of the shift and 
dip model have been compared to another holographic 
model as a baseline comparison. This has illustrated 
that the added phcnomcnological constraints imposed in 
constructing the shift and dip model simultaneously im- 
proves the agreement with the QCD sum rules. Clearly 
this agreement is enhanced by the fact that the moments 
have a dominant contribution from the J / ip and ip' states. 
However, the agreement, in particular for the first mo- 
ment when the contributions from J '/ip and %p' to the 
moment are not as significant, is nontrivial. 

As one progresses to future holographic models of 
heavy-quark systems, it will be important to understand 
how these systems scale with the quark mass. In this way, 
one may be able to relate aspects of the holographic mod- 
els found in heavy-quark systems back to their light mass 
cousins. Moreover, most of the discrepancies between the 
shift and dip model and the QCD sum rules seem to stem 
from discrepancies associated with the first moment. Any 
future attempts to improve the agreement must address 
this issue. Since each moment from the shift and dip 
model is smaller than from the QCD OPE, to improve 
agreement one must either increase the decay constants 
of the excited states or decrease their masses. As dis- 
cussed in the construction of the shift and dip model in 
[l|, these adjustments may be accomplished by changing 
the holographic potential, Eq. (IT51) . by further increasing 
the region of attraction, both in strength and width, in 
the area of z ~ m~ , or by possibly softening the soft 
wall, respectively. Obviously, any future changes would 
still need to be consistent with the spectroscopic data of 
the first two charmonium states. 

In conclusion, we have demonstrated that the holo- 
graphic shift and dip model is consistent in a nontrivial 



7 



manner with the QCD sum rules at zero temperature. 
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Appendix A: Rescaled p model 

For the rescaled p model, the moments can be calcu- 
lated using the procedure outlined at the end of Sec. ITT1 
However, one does not need to solve any equations for 
the function B(z) since it is provided from the soft-wall 
model. This function can be expressed as 



e B(*) = i e -aV 



(Al) 



where a is a parameter of the model which is determined 
by requiring the model reproduce the mass of J /ip cor- 
rectly. Numerically it is a = 1.56GeV. From this ex- 
pression for B(z), one can determine V(z) subject to the 
boundary conditions as 



V(z,q) = T(l 



4a 2 



U 



4a 2 



,0,aV 



(A2) 



where U (a, b, x) is the Tricomi confluent hypergeometric 
function. This expression for V(z,q) can then be used in 
Eq. (fT2j) to determine the polarization function for this 
model. This results in 



n (c) (Q 2 ) = -A 



2 5 I 



(A3) 



where Q 2 = —q 2 , 7e is the Euler gamma, and tp(x) is the 
polygamma (or digamma) function. This function has 
already been renormalized so n( c )(0) = 0. The moments 
of the polarization function can be easily found using 
Eq. Q resulting in, 



1 



1 



2g 2 (4a 2 )« 
where (,(x) is the Zeta function. 



C(n + 1), 



(A4) 



Appendix B: Determining the moments of the 
polarization function for a holographic model 



this may not be the case. In this Appendix, we will cal- 
culate the moments for this arbitrary holographic model 
in terms of integral recursive relations. These can be used 
when an analytic solution is possible, but will most likely 
be more useful when one is required to rely on numerics. 

To construct the moments, we will once again use the 
same procedure that was outlined at the end of Sec. [TT1 
However, since we wish to consider a large class of theo- 
ries, we will leave the function B(z) unspecified. There- 
fore, we begin with the equation of motion for the vector 
field of Eq. 0, 

V"(z, Q 2 ) + B'(z)V'(z, Q 2 ) - Q 2 V(z, Q 2 ) = 0, (Bl) 

where q 2 = —Q 2 . As before, this equation is subject 
to the boundary conditions of Eq. (ITT|) . Since we are 
interested in the polarization function at small Q , we 



can explore Eq. (|B1[) in the same limit. Therefore we 
will assume that V(z, Q 2 ) can be expanded in a Taylor 
series of Q 2 , namely, 



F(z,Q 2 ) = ^(-Q 2 ) I V n (z). 



(B2) 



n=0 



The functions V n (z) correspond to the coefficients in the 
Taylor series. Note these should not be confused with 
the eigenmodes of the differential equation. By inserting 
Eq. (|F32|) into Eq. (IB1[) and equating powers of Q 2 , one 
arrives at a difference equation for the functions V n (z): 



V£ + BX = -Vn-l, 

which can be solved by iterative means. 
The 0th order expression, 



can be easily solved, 

Vo(z) = h 



B% = 0, 



e- B ^dz'. 



(B3) 



(B4) 



(B5) 



The integration constants, b\ and 62, can be determined 
from the boundary conditions resulting in Vq(z) = 1 . 
For larger values of n, the equation can be solved by a 
Green's function technique. Therefore the solutions can 
be written as 



V n (z) 



dz'G(z,z')K-i(A 



(B6) 



where G(z, z') is defined to satisfy the differential equa- 
tion, 

d 2 G(z,z>) , dG(z,z') 

+ B(z) — =-5{z-z). (B7) 



dz 2 



In Sec. IIII1 the moments of the polarization function 
were calculated for a particular holographic model of 
charmonium. In each case, one could solve the necessary 
differential equation analytically, so an analytic expres- 
sion for the polarization function was possible. However, 
for an arbitrary holographic model with a vector field 



3 Though this expression does not explicitly satisfy the boundary 
condition for V(z,Q 2 ) at infinity, we have assumed that V n (z) 
should at least be finite in this limit. The boundary condition for 
V(z, Q 2 ) would then be satisfied once the infinite sum in Eq. {B2j 
is performed. 
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To be consistent with the boundary conditions on 
V(z, Q 2 ) and the expression for Vq(z), G(z, z') must sat- 
isfy the boundary conditions G(e, z') = and G(z, z') is 
hnitc as z — > oo. 

Equation (|F37|) can be solved to the left and to the right 
of delta function as 

I 05 + 06 J z , e { 'dz z > z 

The integration constants can be determined from the 
boundary conditions and the matching conditions at the 
delta function resulting in the Green's function to be ex- 
pressed as, 

G{z, z') = e B ^ [ Z< e^ 2 " W, (B9) 



where z < is the smaller of z and z' . It is useful to note 
the derivative of the Green's function with respect to z, 



dG(z,z') 
dz 



(BIO) 



where 9 is the unit step function. 

Having determined the appropriate Green's function, 
V n (z) for n > 1 can be generated from Eq. (|B6|) . in par- 
ticular V\ (z) can be written as, 



Vi = I G{z,z')V (z')dz' 

' ' e Bi - z ") J' e- B ( z "Uz'^J dz' 
+ J^e B ^ J\' B ^"Uz") ,/:' 

; roc 

dz'e- B ^ / e B ^ z "Uz". 



(Bll) 



The notation here can be compactihcd some by defining 
the functions a(z) and /3(z) as 



a(z) = J e- B(z 'Uz', 
P{z)= I e Bi - z,) dz'. 

J Z 

Using this notation, Eq. (|B9[) can be written as 



(B12) 



= J -I3'(z')a(z) z<z' 
1 ' ' \ -(3'{z')a(z r ) z> z' 



and Vi (z) as 



Vi(z) = I p{z')a'{z')dz' 



(B13) 



(B14) 



Having determined the functions V n in terms of the 
Green's function and V n -i, let us turn our attention to 
the polarization function and its moments. As before, 
Eq. (TT2"j) relates the polarization function to the vector 
field. Using the expansion in Q 2 for V(z, Q 2 ) in Eq. (|B2[) . 
we can expand n( c ) in powers of Q 2 , 



n (c) 4E (-QT (e B(z K + i(z) 



!>:, 



(B15) 



One may worry since this expansion seems to have a 
(Q 2 ) -1 term. However, upon further inspection, the co- 
efficient of this term is b b ^ z 'Vq(z)\ , and with the pre- 
vious assessment that Vq[z) = 1, this term is exactly 
0. The expansion also contains a term constant in Q 2 . 
This term can always be eliminated by renormalizing the 
polarization function such that n( c )(0) = 0. 

The moments can be calculated in the normal manner 
from Eq. ((4]) and yield 



.M n = L( e B {Z ) v , +i{z) \ 

(7c; \ / ; 



95 

% 
1 

~7 5 



G'{e,z')V n (z')dz' 
e B ^V n (z')dz', 



(B16) 



where Eqs. (|B6[) and (|B10|) are used for lines two and 
three respectively. The moment can be expressed using 
the compact notation as 



1 r°° 

JZ n = — [3{z)V;(z)dz. 
95 Jt 

The first moment can be explicitly calculated as 



(B17) 



1 

9% u e 

1 /' ,^-> 



e B ^V 1 (z)dz 



9l „, 



9% 



-B(z') e B(*") dz ») dz > 



dz 



dze' B ^ 



or in the compact notation, 



- / P(z) 2 a'(z)dz. 

91 Je 



(B18) 



(B19) 



All higher moments can then be found by iteratively solv- 
ing Eq. (jB3| for V n (z) and using it in Eq. (|B16j) . One 
can easily see that, for this class of holographic models, 
is positive for all values of n. 
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